ORIGINAL ARTICLE 

http://dx.doi.org/1 0.4306/pi.2014.1 1 .2. 1 92 



Print ISSN 1738-3684 / On-line ISSN 1976-3026 

OPEN ACCESS 



Genetic Role of BDNF Val66Met and 5-HTTLPR Polymorphisms 
on Depressive Disorder 

Kyu Young Lee^ ^ Seong Hoon Jeong^'^, Se Hyun Kim", Yong Min Ahn"'^'^'', Yong Sik Kim^, 
Hee Yeon Jung^^ Yang Weon Bang^", and Eun-Jeong Joo^'^^ 

^Department of Neuropsychiatry, Eulji University School of Medicine, Daejeon, Republic of Korea 
^Department of Psychiatry, Eulji General Hospital, Seoul, Republic of Korea 
'Department of Psychiatry, Eulji University Hospital, Daejeon, Republic of Korea 

"Institute of Human Behavioral Medicine, Medical Research Institute, Seoul National University, Seoul, Republic of Korea 
^Department of Psychiatry and Behavioral Science, Seoul National University College of Medicine, Seoul, Republic of Korea 
''Department of Neuropsychiatry, Seoul National University Hospital, Seoul, Republic of Korea 
'Biomedical Research Institute, Seoul National University Hospital, Seoul, Republic of Korea 

''Department of Neuropsychiatry, Dongguk University International Hospital, Dongguk University Medical School, Goyang, Republic of Korea 
"Department of Psychiatry, SMG-SNU Boramae Medical Center, Seoul, Republic of Korea 
^"Department of Psychiatry, Keyo Hospital, Uiwang, Republic of Korea 



Objective We investigated possible association between depressive disorders and BDNF Val66Met and 5-HTTLPR. Brain derived 
neurotrophic factor (BDNF) gene and serotonin transporter (SLC6A4) gene are promising candidate genes for depressive disorders. It 
has been suggested that BDNF promotes the survival and differentiation of serotonergic neurons and that serotonergic transmission ex- 
erts powerful control over BDNF gene expression. 

Methods Final analyses were performed on 186 patients with depressive disorders and 1032 controls. Val66Met polymorphism of 
BDNF gene and 5-HTTLPR polymorphism of serotonin transporter gene were genotyped and allele and genotypic associations on the 
diagnosis of depression and age at onset of depression were analyzed. 

Results The 5-HTTLPR was positively associated with depressive affected status in the total sample and in females (p=0.038 for allele- 
wise, p=0.015 for genotype-wise associations), but, not in males. The BDNF Val66Met showed no association with depression. BDNF 
Val66Met and 5-HTTLPR alone were not associated with age at onset of depression. Additional analysis on the interaction between 
BDNF Val66Met and 5-HTTLPR found a significant association with age at onset of depression in the entire patient group. This associa- 
tion was also found in the female but not in the male patient group. None of the positive results survived Bonferroni correction for mul- 
tiple testing. 

Conclusion This result suggested that BDNF Val66Met and 5-HTTLPR may contribute to depressive disorders in a complex way and 
that the genetic effect could differ by gender. Further studies with large number of patients will be necessary. 
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INTRODUCTION 

Depressive disorder is a complex disease phenotype. Genet- 
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ic Studies have indicated that genes work significantly as etio- 
logical factors and disease modifiers in depression. However, 
no major genes for depressive disorder likely exist. Accumulat- 
ing genetic studies such as candidate gene studies, genome- 
wide association studies, and rare-variant studies strongly sug- 
gest that polygenes with small effects and rare mutations are 
likely genetic precipitating factors. Among various promising 
candidate genes for depression, the serotonin transporter (SL- 
C6A4) and brain-derived neurotrophic factor (BDNF) genes 
are two of the most studied candidate genes for depressive 
disorders worldwide. 
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The serotonin transporter-linked polymorphic region 
(5HTTLPR) for SLC6A4 is located upstream of the transcrip- 
tion initiation site of the SLC6A4 gene. 5HTTLPR is a 44-bp 
insertion/deletion variable-number tandem repeat polymor- 
phism with two common alleles: the 16-repeats long allele 
and the 14-repeat short allele. 5HTTLPR is a functional poly- 
morphism in that the short allele has lower transcriptional ef- 
ficiency.''^ Studies on 5-HTTLPR and depressive disorders in- 
cluding meta-analyses can be summarized as follows. First, 
the association with depression itself seems to be positive. 
Clarke et al.^ reported a small but statistically significant asso- 
ciation based on a meta-analysis of 39 studies. Second, an as- 
sociation with the antidepressant response is less likely. Taylor 
et al. analyzed 28 studies and found no significant effect on 
remission rates after correcting for missing data.'' Third, the 
genetic moderation by 5-HTTLPR of environmental factors 
for the development of depression is strongly supported. Ear- 
lier studies found a significant interaction between 5-HTTL- 
PR and the number of life depressive events, and the short al- 
lele indicates a higher vulnerability to social stressors. This 
suggests that environmental insult is moderated by this gene.^'*" 
A meta-analysis including 14 studies found no evidence of an 
association between the serotonin-transporter genotype alone 
and no interaction with stressful life events.' However, Karg et 
al.* recently included 54 studies in a meta-analysis and found 
that the s allele was associated with an increased risk of de- 
pression under stress. 

BDNF is one of the neurotrophic factors manifested in the 
brain and multiple tissues and has been suggested as a bio- 
marker of depressive disorders.' Findings indicating that 
BDNF levels are reduced in the postmortem brain and blood 
of patients with depressive disorders and that antidepressant 
treatments increase BDNF levels support the neurotrophic 
hypothesis for the etiology of depression.'"" BDNF Val66Met 
(G196A, rs6265) polymorphism is located at an exon and 
changes an amino acid at codon 66 from Val to Met. Positive 
association studies have been published,'^ ''' however, whether 
the Met or Val allele is a risk allele for depression has not been 
confirmed. FrieUngsdorf et al. reported that Met homozygotes 
may have a greater risk for major depressive disorder ''' An- 
other study found that Val homozygotes have an increased 
chance of depression." However, negative association results 
have also been found,'^"" and a meta-analysis involving 14 
studies reported that BDNF Val66Met was not significantly 
associated with depression in a total sample and that a signif- 
icant effect was found only in men." Environmental interac- 
tions with this SNP that increase susceptibility for major de- 
pressive disorder after early-life stress may be associated with 
the Met allele.'* Taken together, these results indicate that the 
genetic effect of BDNF Val66Met on depression is not fuUy 



understood. 

A significant interaction between BDNF and serotonin has 
been suggested in brain function. Based on a common feature 
of their regulatory ability for the development and plasticity 
of neural circuits, the interaction of these two systems is like- 
ly involved in mood disorders. BDNF promotes the survival 
and differentiation of serotonergic neurons, and serotonergic 
transmission exerts powerful control over BDNF gene ex- 
pression." Administration of SSRIs positively regulates BDNF 
transcription through increases in CREB phosphorylation.^" 
Other evidence for this positive interaction is substantial. 
However, specific combinations of risk alleles are not consis- 
tent. The Met allele of the BDNF Val66Met polymorphism 
protects against deleterious effects of the 5-HTTLPR s allele 
on brain volume."' The Met allele seems to confer a better re- 
sponse to lithium in s-carrier patients.^^ In contrast, the BDNF 
Met allele carriers and 5-HTTLPRs homozygotes have the 
highest depression vulnerability under adverse environ- 
ments.^^'^"* Understanding the interactions between these two 
systems is crucial to help obtain a more comprehensive view 
of the molecular pathways of mood regulation." 

Age at onset is an important clinical feature of disease and 
may be influenced by genetic factors. Genetic factors for the 
age at onset could be identical to or different from the genetic 
factors in the diagnosis of depressive disorder In the present 
study, we investigated a possible association between depres- 
sive disorder and BDNF Val66Met and 5-HTTLPR, which 
are promising candidate genetic polymorphisms in psychiat- 
ric genetics. The influences of individual genes and of the in- 
teraction between polymorphisms on the diagnosis of depres- 
sion and age at onset were explored. The association between 
age at onset of depressive disorder and a genetic interaction 
between BDNF Val66Met and 5-HTTLPR has not been in- 
vestigated previously. This study aimed to replicate the genetic 
association of BDNF Val66Met and 5-HTTLPR with depres- 
sive disorder in a Korean population and to explore the inter- 
active genetic role of BDNF Val66Met and 5-HTTLPR on age 
at onset of depressive disorders. 

METHODS 
Subjects 

Patients were recruited from Eulji General Hospital. All pa- 
tients were evaluated based on the diagnostic criteria of the 
DSM-IV. We included a broad spectrum of depressive disor- 
ders for this study, i.e., major depressive disorder, dysthymia, 
and depressive disorder, not otherwise specified (NOS). Both 
single-episode and recurrent major depressive disorder were 
included, and major depressive disorder with dysthymia was 
also included. Patients were individually interviewed by re- 
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search nurses using the Korean version of the Mini-Interna- 
tional Neuropsychiatric Interview. Consensus diagnostic 
meetings of more than two psychiatrists, including the psy- 
chiatrist in charge of the psychiatric treatment of the patient, 
made the final diagnosis of participants. Associated symp- 
toms and other clinical characteristics (Table 1) were collected 
based on patient-provided information. Psychotic features, 
psychomotor retardation, somatization, suicide attempts, 
panic symptoms, alcohol problems, and familial loading of 
major psychiatric disorders were evaluated by the psychiatrist 
in charge. Evaluation of weight gain, increased appetite, and 
increased sleep were largely based on subjective information 
from patients. Subjects with a history of any kind of organic 
abnormality of the brain, substance dependence, drug abuse, 
or other physical conditions possibly manifesting as psychi- 
atric symptoms were excluded from this study. Controls were 
recruited from college students, nurses, and fire and public 
protection officers. A brief psychiatric interview was per- 
formed by psychiatric research nurses to evaluate current and 
past psychiatric illness. Subjects with a lifetime history of ma- 
jor psychiatric illness and/ or brain trauma were excluded. Fi- 
nal analyses were performed on 186 patients and 1032 con- 
trols. The patient group consisted of 106 patients with re- 
current major depressive disorder, 38 patients with a single 



episode of major depressive disorder, 17 patients of major de- 
pressive disorder with dysthymia, seven patients with dysthy- 
mia, and 19 patients with depressive disorder, NOS. Control 
subjects, matched for sex ratio with the patient group, were 
randomly selected from our control pool. The mean ages of 
the patient and control groups were 52.14±14.38 and 24.43± 
5.07 years, respectively. More detailed information regarding 
the subject number and mean age for each group is given in 
Table 2. All subjects participating in this study signed a writ- 
ten informed consent form, and the study protocol was ap- 
proved by the ethics committee of Eulji General Hospital. 

Genotyping 

DNA was extracted from blood samples using a DNA iso- 
lation kit (Roche, Mannheim, Germany). The Val66Met SNP 
of the BDNF gene and the 5-HTTLPR polymorphism of the 
serotonin transporter gene were genotyped. Genotyping for 
SNP Val66Met was performed using the TaqManTM method 
(Applied Biosystems, Foster City, CA, USA).^^ Primers for the 
Val66Met G/A polymorphism were synthesized with the se- 
quences 5'-TCAAGAGGCTTGACATCATTGG-3' for the 
forward primer and 5'-GCCGAACTTTCTGGTCCTCAT-3' 
for the reverse primer FAM and VIC dye probes were also 
synthesized with the sequences 5'-6FAM-ACACTTTCGAA 



Table 1 . Clinical characteristics of patients with depression 



Total 
Number 



Male 



Number 



% 



Female 
Number % 



Diagnosis 
MDD signle 
MDD recurrent 
MDD with dysthymia 
Dysthymia 

Depressive disorder, NOS 
Associated symptom 

Psychotic feature 

Psychomotor retardation 

Weight gain 

Increased appetite 

Increased sleep amount 

Somatization 

Suicidal attempt 

Panic symptom 
Characteristics 

Alcohol problem 

Familial loading 

Postmenopausal status 



37 
106 
17 
7 
19 

23 
56 
26 
14 
11 
151 
22 
26 

8 
46 



11 
19 
3 
3 
3 

10 
14 
3 
1 
2 
25 
6 
1 

6 
9 



28.21 
48.72 
7.69 
7.69 
7.69 

25.64 
35.90 
7.69 
2.56 
5.13 
64.10 
15.38 
2.56 

15.38 
23.08 



26 
87 
14 
4 
16 

13 
42 
23 
13 
9 
126 
16 
25 

2 
37 
87 



17.69 
59.18 
9.52 
2.72 
10.88 

8.84 
28.57 
15.65 
8.84 
6.12 
85.71 
10.88 
17.01 

1.36 

25.17 
59.18 



MDD: major depressive disorder, NOS: not otherwise specified 
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CACGTGATAGAAGAGCTGTT-TAMRA-3' for the FAM 
dye probe and 5'-TET-ACACTTTCGAACACATGATAG 
AAGAGCTGTTG-TAMRA-3' for the VIC dye probe. Stan- 
dard polymerase chain reaction (PGR) analyses was carried 
out in 5-\iL volumes that contained 50 ng of genomic DNA, 
2.5 i^L of 2x TaqMan Master Mix, 600 nM of each primer, 
and 100 nM of FAM and VIC dye probes. After an initial 2 
minutes of denaturation at 50°C followed by 10 minutes at 
95°C, 40 thermal cycles consisting of 15 seconds at 95°C and 
1 minute at 60°C were carried out. Automated reading was 
done with an ABI Prism 7900 HT Sequence Detection System 
(Applied Biosystems). The 5-HTTLPR polymorphism was 
genotyped as follows. Primers for 5-HTTLPR were synthe- 
sized with the sequences for stpr5: 5'-GGCGTTGCC 
GCTCTGAATGC-3' and stprS: 5'-GAGGGACTGAGCTG 
GACAACCAC-S: A l-|iL aliquot of 50 ng/^L DNA, 2 ^iL of 
20 mM dNTP mix (Roche), 0.5 [iL of each 20 pM primer, 1.5 
[jL of lOx buffer, 1.5 U of Taq polymerase (Roche), and 9.2 
\iL of distilled water were mixed to make a 15-[jL PGR reac- 
tion volume. The PGR was carried out using the following 
steps: 5 min at 94°G to denature the DNA sequence, 35 cycles 
at 94°G for 30 sec, 64°G for 30 sec, and 72°G for 1 min, fol- 
lowed by 7 min for amplification at 72°G. PGR products were 
examined by electrophoresis for 30 min on a 3% agarose gel 
stained with ethidium bromide. DNA sequencing of the PGR 
products for several subjects was performed to confirm the 
5-HTTLPR genotype. 

Statistical analysis 

Hardy-Weinberg equilibrium was tested by a goodness-of- 
fit chi-square test. The entire analysis was conducted using 
UNPHASED 3.0.13 software (http://www.mrc-bsu.cam.ac. 
uk/personal/frank/software/unphased/). Allelic and genotyp- 
ic associations of each polymorphism were tested for the as- 
sociation with diagnosis of depression. The program calculat- 
ed chi-square values and p- values for allelic and genotypic 
interactions according to a full model. The level of statistical 
significance was set at p=0.05. Permutation tests were per- 
formed 1000 times to avoid false discovery from initial analy- 
ses. Age at onset was tested for an association with each poly- 
morphism by the chi-square test with UNPHASED 3.0.13 



and analysis of variance (ANOVA). Three groups for BDNF 
Val66Met (Met/Met, Met/Val, and Val/Val) and two groups 
for 5-HTTLPR (ss and si, and 11) were constructed for six 
combined-genotype groups. Analysis for a possible associa- 
tion between age at onset and an interaction of BDNF Val- 
66Met with 5-HTTLPR was performed by ANOVA for these 
six groups. The ANOVA was performed using SPSS 12.0 soft- 
ware (Chicago, IL, USA). 

RESULTS 

The patient group consisted of 147 females and 39 males. 
As females were predominant in the patient group, we ran- 
domly selected control subjects under the condition of match- 
ing the sex ratio with the patient group. Major depressive dis- 
order was the most common diagnosis in this patient group. 
Major depressive disorders constituted 77% of the entire pa- 
tient group. About 60% of patients were recurrent major de- 
pression and their age and age at onset were not significantly 
different from the other patients. The age at onset was 44.02± 
15.32 years old for the total patient group and was not statisti- 
cally different according to gender (Table 2). Somatization 
symptoms and panic symptoms associated with depressive 
disorder were found significantly more often in the female 
patient group (p=0.005 for somatization and p=0.019 for panic 
symptoms; chi-square test). Associated psychotic features and 
alcohol problems were significantly more frequent in the male 
patient group (p=0.011 for psychotic features and p<0.001 for 
alcohol problems, chi-square test). More than half of the fe- 
male patients (59%) were postmenopausal (Table 1). 

The serotonin transporter 5-HTTLPR polymorphism was 
positively associated with depression. Both allele and geno- 
type associations were significant in the total sample (p=0.038 
for allele-wise, p=0.015 for genotype-wise associations). 
When we conducted the analysis separately according to gen- 
der, a significant genetic association was found in the female 
group but not in the male group. The significance survived 
1000 permutations for both the total and female groups (Table 
3). The BDNF Val66Met SNP was not associated with depres- 
sion (Table 4). 

Age at onset as a quantitative phenotype was tested for a 



Table 2. Age and age at onset of subjects 







Depression 


Control 




Recurrent depression 




Total 


Male 


Female Total Male 


Female 


Total Male Female 


Number 


186 


39 


147 1032 219 


813 


106 19 87 


Age 


52.14(14.38) 


52.67 (17.73) 


52.00 (13.42) 24.43 (5.07) 27.10 (8.15) 


23.71 (3.52) 


54.87(11.90) 54.74(17.36) 54.90(10.48) 


Age at onset 


44.02 (15.32) 


45.85 (19.62) 


43.53 (13.99) 




44.56(13.41) 44.58(18.80) 44.55(12.07) 



Mean values for age and age at onset were presented. SD were pressented in parenthesis 
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possible association with BDNF Val66Met, 5-HTTLPR, and 
the interaction of polymorphisms. Both the chi-square and 
ANOVA tests consistently found that neither BDNF Val- 
66Met alone nor 5-HTTLPR alone was associated with age at 
onset of depressive disorder in our sample (Table 3 and 4). The 
Met/Met genotype of the BDNF Val66Met polymorphism was 
associated with the oldest age at onset (46.65± 16.61 years), 
and the Val/Val genotype with the youngest age at onset 
(41.81±14.98 years). The Met allele tended to be associated 
with older age at onset in a dose-dependent manner, and this 
association was significant (p=0.293) (Table 5). Both the ss 
genotype and si genotype of the 5-HTTLPR polymorphism 
showed similar age at onset, and these tended to be different 
from the 11 genotype. The association between the 5-HTTLPR 
genotype and age at onset did not differ significantly among 
the three genotypes (p=0.705) (Table 5). We found a signifi- 
cant effect of the interaction between BDNF Val66Met and 
5-HTTLPR and the age at onset of depression (Table 6). The 
positive result in the total patient group (F=2.808, df=5, 
p=0.018) (Table 6), and in the female group (F=2.379, df=5, 
p=0.042, data not shown), but not in the male patient group 
(F=1.413, df=5, p=0.245, data not shown) were revealed. 

DISCUSSION 

As in many studies, a genetic association between 5-HT- 
TLPR and depressive disorder was found in this Korean sam- 
ple. The s allele frequency in control group was 77.52%, simi- 

Table 3. Genetic association between BDNF Val66l\/let and depression 



lar to other studies on Korean population.^^'" We found a 
deviation from the Hardy- Weinberg equilibrium for 5-HT- 
TLPR in the total and female patient subgroup but not in the 
male subgroup (Table 3). This result could add more support 
for the positive association between 5-HTTLPR and depres- 
sive disorder, along with the chi-square test for the allele and 
genotype association, although many factors influence Hardy- 
Weinberg equilibrium. This association was not found in the 
male subgroup. It is probably because of the small number of 
subjects. We included a relatively large number of control 
subjects for the comparison to overcome the statistical weak- 
ness based on the small number of patients in this study. An- 
other reason could be a gender difference in genetic modula- 
tion of depression by 5-HTTLPR. A gender-specific genetic 
effect has been reported for 5-HTTLPR in many psychiatric 
phenotypes, particularly in terms of the interaction with an 
environmental effect.^*'^' It is possible that the genetic effect is 
moderated in a gender-specific way. However, considering 
the limited statistical power of this study, further studies with 
a larger sample size are necessary. 

BDNF Val66Met did not show an association with depres- 
sive disorder, and the interaction between BDNF Val66Met 
and 5-HTTLPR was found not to be associated with the diag- 
nosis of depressive disorder in this study. The A allele (Met) 
frequency was 46.66% in this control group, not significantly 
different from other studies on Korean population.^'''^^'^" There 
was no deviation from Hardy- Weinberg equilibrium for 
BDNF Val66Met in our sample (Table 4). Despite much bio- 



Total 



Male 



Female 



Allele 
1 2 



Genotype 
11 12 22 



Allele 
1 2 



Genotype 
11 12 22 



Allele 
1 2 



Genotype 
11 12 22 



Frequency 

Depression 

Control 
HWE 

Depression 

Control 
Diagnosis of depression 

p-value 

Permutation p-value 
Age at onset 

p-value 

Permutation p-value 



172 200 
963 1101 



43 
218 



86 
527 



57 
287 



40 
212 



38 
226 



10 
54 



20 
104 



9 
61 



X^=0.91 1 , p-value=0.634 
X^=0.692, p-value=0.707 



X^=0.027, p-value=0.987 
X'=0.531,p-value=0.767 



132 162 
751 875 

v2— 



33 66 
164 423 



48 
226 



X^= 1.260, p-value= 
X^=1.771, p-value= 



=0.533 
=0.412 



0.022 
0.881 
0.863 

2.657 
0.103 
0.108 



1.472 
0.479 
0.489 

2.493 
0.286 
0.281 



0.220 
0.639 
0.642 

0.926 
0.336 
0.349 



0.401 
0.818 
0.826 

3.530 
0.171 
0.191 



0.167 
0.683 
0.654 

1.549 
0.213 
0.247 



2.583 
0.275 
0.268 

2.792 
0.248 
0.274 



HWE: Hardy- Weinberg Equilibrium, Permutation p-value: p-value of permutation test in UNPHASED 3.0.13 (1000 times), BDNF: brain de- 
rived neurotrophic factor 
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Table 4. Genetic association between 5-HTTLPR and depression 




Total 






Male 




Female 






Allele Genotype 




AUele 


Genotype 


Allele Genotype 




1 2 11 


12 22 




1 2 11 


12 22 


1 2 11 


12 22 


Frequency 
















Depression 


306 66 133 


40 13 




64 14 28 


8 3 


242 52 105 


32 10 


Control 


1600 464 636 


328 68 




357 81 149 


59 11 


1243 383 487 


269 57 


HWE 
















Depression 


X'=12.887, p-value=0.002 




X'=3.594, p-value=0.166 


X^=9.364, p-value=0.009 


Control 


X^=8.011,p-value=0.018 




X'=2.476, p-value=0.290 


X^=5.366, p-value=0.068 


Diagnosis of depression 


















4.327 


8.413 




0.013 


1.022 


5.130 


8.177 


p-value 


0.038 


0.015 




0.909 


0.600 


0.024 


0.017 


Permutation p-value 


0.050 


0.015 




0.884 


0.592 


0.048 


0.015 


Age at onset 


















0.267 


0.714 




0.136 


0.914 


0.823 


0.738 


p-value 


0.606 


0.700 




0.712 


0.633 


0.364 


0.691 


Permutation p-value 


0.658 


0.700 




0.766 


0.630 


0.409 


0.691 


HWE: Hardy-Weinberg Equilibrium, Permutation p-value: p-value of permutation test 


in UNPHASED 3.0.13 (1000 times) 




Table 5. Age at onset for genotypes of BDNF Val66IVIet and 5-HTTLPR 








Gene 


Genotype 




N 


Mean 


SD 


F 


p-value 


BDNF VaI66Met 


Met/Met 




43 


46.65 


16.61 


1.236 


0.293 




Met/Val 




86 


44.16 


14.82 








Val/Val 




57 


41.81 


14.98 






5-HTTLPR 


ss 


133 


44.13 


15.84 


0.350 


0.705 




si 




40 


44.73 


14.13 








11 




13 


40.69 


13.88 






ANOVA test was performed. BDNF: brain derived 


neurotrophic factor, 5-HTTLRP: The serotonin transporter-linked polymorphic region 


Table 6. Genetic interaction of 5-HTTLPR and BDNF Val66Met on the age at onset of depression 
















5-HTTLPR 














ss & si 




U 








Number 




Mean 


SD Number Mean 


SD 


BDNF VaI66Met 


AA (Met/Met) 


37 




49.11 16.39 


6 31.50 


7.97 




AG (Met/Val) 


83 




44.36 15.04 


3 38.67 


1.16 




GG (Val/Val) 


53 




40.74 14.66 


4 56.00 


13.37 



ANOVA test was performed (F=2.808, df=5, p=0.018). 5-HTTLPR serotonin transporter-linked polymorphic region, s: short allele, 1: long 
allele, BDNF: brain-derived neurotrophic factor. A: alanine, G: guanine. Met: methionine, Val: valine, SD: standard deviation 



chemical evidence for the role of BDNF in depressive disor- 
ders from studies such as animal models of depression, blood 
samples, and brain tissue of patients with depression, previous 
genetic association studies have not been consistent with re- 
spect to the association between BDNF Val66Met and depres- 
sive disorder. Even for Korean population, positive associa- 
tion with poststroke depression and depression in women 
with mastectomy due to breast cancer^'"'^' as well as no associ- 



ation with major depression™ have been reported. Our find- 
ings do not support a conclusion that the BDNF gene plays a 
significant role in depressive disorder in Korean population. 

A gene-gene interaction between 5-HTTLPR and the 
BDNF Val66Met polymorphism was found on the age at on- 
set of depressive disorder. It is unclear which allele conferred 
a greater risk for depressive disorder in this study. The Met al- 
lele of BDNF Val66Met was associated with older age at onset 
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of depressive disorder if the patient had either the 5-HTTLPR 
ss or si genotype. However, the BDNF Val66Met Val allele was 
associated with older age at onset of depressive disorder in the 
patient group with the 11 5-HTTLPR genotype, indicating that 
the BDNF Val66Met polymorphism affects age at onset of de- 
pressive disorders insufficiently by itself but instead interacts 
with 5-HTTLPR. Early onset depressive disorder is consid- 
ered to be different from late onset depression in some as- 
pects. Genetic attribution to the age at onset is most likely to 
exist and several genes, not only one gene may play roles in 
interaction. 

Genetic epistasis between BDNF Val66Met and 5-HTTLPR 
has been investigated for various psychiatric phenotypes. 
5-HTTLPR 11 carriers scored lower on neuroticism scales in 
the presence of the BDNF Val variant and higher in the pres- 
ence of the BDNF Met variant.'' The 5-HTTLPR ss genotype 
had a negative impact on mental health after childhood abuse 
only in patients with the BDNF Val/Val genotype.*" A gene- 
gene interaction on conscientiousness in adolescence was 
found in BDNF Val66Met Met allele carriers, with partici- 
pants with the 5-HTTLPR ss genotype having the lowest con- 
scientiousness scores.'^ These studies support that a 5-HTTL- 
PR and BDNF Val66Met gene-gene interaction plays a role 
on the phenotypes related to mood disorders and personality. 
However, no specific combination of risk alleles has been 
consistently found. Thus, replication studies with a larger 
sample size are warranted. A positive influence of the interac- 
tion between 5-HTTLPR and BDNF Val66Met on age at on- 
set of depressive disorder was found in the total sample and 
in the female patients but not in the male patients. It is possi- 
ble that we could not detect the difference statistically because 
our male subgroup was too small. Another reason for this 
could be a gender difference in the nature of the genetic asso- 
ciation. Many studies have reported a gender difference in the 
genetic association of BDNF Val66Met''''' and 5-HTTL- 
PR"' with depression. Clinically, depressions are clearly dif- 
ferent according to gender in many ways. More studies focused 
on gender specific gene-gene interaction will be necessary. 

This study had the following limitations. The number of 
patients included in the study was insufficient and far fewer 
than the number of control subjects. This led to a statistical 
limitation in interpreting the results in greater depth. Addi- 
tionally, the polymorphisms tested represented only one site 
for each gene: 5-HTTLPR and BDNF. Even though these are 
functional polymorphisms, it is difficult to say whether they 
represent the entire gene. Another statistical limitation came 
from testing of multiple hypotheses. We conducted analyses 
for the whole group and subgroups by gender as well as an al- 
lele and genotype comparison between patients and controls. 
None of the positive results survived the Bonferroni correc- 



tion for multiple testing. Interpretation of these results should 
take into account the fact that the p-values presented in this 
study were not corrected for multiple testing even though the 
permutation tests were performed and our positive results 
still showed significant p-values on permutation tests. Genetic 
power calculation of our sample was done for each polymor- 
phism. Considering the minor allele 1 as a risk allele for 5-HT- 
TLPR and Met allele for BDNF Val66Met and a significance 
level of p=0.05, our analysis had a power of 79.42% for 
5-HTTLPR and 65.23% for BDNF Val66Met under a domi- 
nant model (genotypic relative risk of AA and Aa=1.5; disease 
prevalence=0.10) calculated using the Genetic power calcula- 
tor (http://pngu.mgh.harvard.edu/~purcell/gpc/cc2.html). 
Population stratification due to possibly different underlying 
genetic backgrounds between patients and controls should be 
considered as one of the causes of false-positive results for this 
type of genetic association study. However, the subjects were 
all ethnically Korean; therefore, it is unlikely that population 
stratification occurred. 

In conclusion, we found that 5-HTTLPR, but not BDNF 
Val66Met, contributed to the diagnosis of depressive disorder. 
An interaction between BDNF Val66Met and 5-HTTLPR af- 
fected the age at onset of depressive disorders. Our findings 
suggest that BDNF Val66Met and 5-HTTLPR are significant- 
ly associated with depressive disorders in a complex way, not 
only in the diagnosis of depressive disorder but also in various 
clinical manifestations such as age at onset. The extent of ge- 
netic contribution of BDNF Val66Met and 5-HTTLPR could 
be different between males and females. Understanding the 
interactions between BDNF and the serotonin transporter is 
crucial to obtain a more comprehensive view of the molecular 
pathways of mood regulation. Our results should be consid- 
ered cautiously because of insufficient power based on a rela- 
tively small number of patients. Further studies on this issue 
should follow. 
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